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Abstract 

We consider the relations between the gauge couplings at the electroweak scale 
and the high scale where unification of the three gauge couplings is expected. 
Threshold corrections are incorporated both at the supersymmetric and at the grand 
unified scale and, where available three-loop running and two-loop decoupling are 
employed. We study the impact of the current experimental uncertainties of the 
coupling constants and the supersymmetric mass spectrum on the prediction of the 
super-heavy masses within the so-called minimal supersymmetric SU(5). As a main 
result of the three-loop analysis we confirm that minimal supersymmetric SU(5) 
cannot be ruled out by the current experimental data on proton decay rates. 
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1 Introduction 



An appealing hint in favour of supersymmetry (SUSY) is the apparent unification of gauge 
couplings at a scale of about 10 16 GeV p~H3]. In particular, the gauge coupling unification 
is predicted, even under the assumption of a minimal particle content of the underlying 
Grand Unified Theory (GUT), the so-called minimal SUSY SU(5) model [HE]. An im- 
portant feature of this model is that the requirement of minimality renders it the most 
predictive model from the currently known candidates for GUTs. However, immediately 
after the formulation of the minimal SUSY SU(5) it has been noticed that within SUSY 
GUTs new dimension-five operators cause a rapid proton decay [6j[7]. This aspect was 
intensively studied over the last thirty years with the extreme conclusions of Refs. [HIE] 
that the minimal SUSY SU(5) model is ruled out by the combined constraints from pro- 
ton decay and gauge coupling unification. However, in the following years several careful 
analyses have shown that the constraint on the coloured Higgs triplet mass from proton 
decay was actually too strong, and that minimal SUSY SU(5) is still a perfectly viable 
theory. The proton decay rate for the dominant channel p — > K + V can be suppressed 
either by sfermion mixing [TU] or by taking into account higher dimensional operators in- 
duced at the Planck scale [TTHT3] . Such operators can also lead to a successful down quark 
and charged lepton Yukawa coupling unification for the first two families, which is not 
possible in the renormalizable version of the minimal SUSY SU(5). Planck-scale- induced 
operators also change the boundary conditions for the gauge couplings at the GUT scale, 
since in their presence the colour octet and the isospin triplet Higgs bosons contained in 
the 24# representation do not necessarily have the same mass anymore, as it is the case 
in the minimal model. This fact allows for significantly weaker constraints from gauge 
coupling unification on the colour triplet Higgs boson than in the renormalizable model. 

Recently, new experimental data for the relevant input parameters [T4T4T6] and substantial 
progress on the theory side [TTH20] became available. This encourages us to reanalyze 
minimal SUSY SU(5) focusing on the aspect of gauge coupling unification. We emphasize 
that in our analysis we restrict ourselves to the renormalizable version of minimal SUSY 
SU(5) although in this way the unification of the first and second generation Yukawa 
couplings can not be achieved. In our analysis we ignore this fact since it is not relevant 
here. The constraints on the mass of the coloured Higgs triplet that are derived can easily 
be translated to the non- renormalizable version of minimal SUSY SU(5) by comparison 
of, e.g., Eqs. (9) and (10) of Ref. [13] and Eq. (1T41) in our paper. In this sense, our results 
represent the "worst case scenario". 

More precisely, we review in this paper the constraints on the mass of the coloured Higgs 
triplet Mh c and the grand unification sea la3 Mq within minimal SUSY SU(5) [H|5] taking 
into account the latest experimental data for the weak scale parameters and the most 
precise theoretical predictions currently available. We adopt the renormalization group 
method in the "bottom-up" approach to predict the values of the two parameters and take 

1 See Section [2] for the exact definition of Mh c and Mq. 
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into account threshold corrections generated by the superpartners of the SM particles as 
well as those due to the super-heavy SUSY-GUT particles. In this way, we can derive the 
two SUSY-GUT parameters from the knowledge of the gauge coupling constants of the SM 
at the electroweak scale and the minimal supersymmetric SM (MSSM) mass spectrum. 
In addition, we implement the perturbativity restrictions for all gauge, Yukawa and Higgs 
self couplings, i.e., we require that they are smaller than one up to the Planck scale. 

Besides the minimal SU(5) we also briefly discuss the phenomenological consequences of 
the Missing Doublet Model (MDM) which has been designed in order to avoid unnatural 
doublet-triplet splitting in the Higgs fields of the 5 and 5 representations. 

The remainder of the paper is organized as follows: In the next Section we introduce 
our framework and describe the tools we have used for our analysis. In particular we 
specify the underlying GUT theory and describe in detail our procedure for the running 
and decoupling. In Section |3] the phenomenological consequences are discussed where we 
study in particular the constraints on the GUT masses of minimal SUSY SU(5). Finally, 
we present our conclusions in Section [4j 



2 Theoretical framework 

2.1 Minimal SU(5) and Missing Doublet Model 

In Section [3] we discuss in detail the restrictions on "minimal SU(5)", however, mention 
also briefly the consequences for the so-called "Missing Doublet Model" . For convenience 
of the reader we introduce in this Subsection some details on these models. 

The superpotential of minimal SU(5) [I] is given by 

W = MiTr(S 2 ) + AiTr(S 3 ) + A 2 HEH + M 2 HH 

+ yfiTtfVihH + ^%9jH , (1) 

where and <pi (i = 1,2,3 is a generation index) are matter multiplets in the 10- 
and 5-dimensional representation of SU(5), respectively, and the field H (H) is real- 
ized in the 5 (5) representation. SU(5) is broken to SU(3)xSU(2)xU(l) if the ad- 
joint Higgs boson E = $ a T a (a = 1, . . . ,24) gets the vacuum expectation value (E) = 
U/(2 v / 30)xdiag(-2,-2,-2,3,3), with V = -4V50Mi/(3A 1 ). Choosing (H) = (H) < V 
and in addition imposing the (tree-level-) fine-tuning condition M 2 = —y/ZX-iV/y/AQ the 
isodoublets in H and H remain massless. Furthermore, one gets the following super-heavy 
mass spectrum: 

M 2 x = ^g 2 V\ M 2 Hc = ^X 2 2 V\ Ml ee M ( 2 8jl) = M ( 2 lj3) = 25M 2 il} = l ^\ 2 V 2 , (2) 

where the indices in round brackets refer to the SU(3) and SU(2) quantum numbers. 
Here Ms denotes the mass of the colour octet part of the adjoint Higgs boson E and 
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Mh c stands for the mass of the colour triplets of H and H, M x is the mass of the gauge 
bosons and g is the gauge coupling. The equality M 2 8 ^ = M 2 X 3 ^ only holds if one neglects 
operators that are suppressed by 1/Mpi as we do here. We emphasize again that taking 
into account such operators can considerably weaken the constraint on the colour triplet 
mass Mh c that is derived later in this paper [TTHT3] . However, the altered constraint can 
easily be derived from our results. 

The effects of the super-heavy particle masses on the MSSM gauge couplings can be 
parametrized with the help of decoupling coefficients. If we interpret the MSSM as the 
low-energy effective theory of the SUSY SU(5) model, we can define its three gauge 
coupling constants as functions of the unique SU(5) gauge coupling a su ( 5 ) through 

af SSM (/i G UT) = CJ/icuT, « SU(5 Uh c , M x , M e ) a su(5) C"gut) , z = 1, 2, 3 , (3) 

where Co^^gut, a su ^ 5 \ M Hc , M x , Ms) denote the decoupling coefficients that can be ob- 
tained from Green's functions with external light particles computed in the full and ef- 
fective theory. The scale /igut is an unphysical parameter, not fixed by theory. The 
dependence of physical observables on /igut thus provides an estimation of the theoreti- 
cal uncertainties within fixed order perturbation theory. 

For SUSY theories the most convenient regularization scheme is Dimensional Reduction 
(DRED) [2 lj which we also adopt in our calculation. As a consequence, the coupling 
constants appearing in Eq. ([3]) are renormalized minimally in the so-called DR renormal- 
ization scheme. The one-loop formulas of the decoupling coefficients for a general gauge 
group have been known for a long time [221425] . The specification to minimal SUSY SU(5) 
reads 



a sv ^(u) ( 2 
C*(A0 = 1+ An W l~L, Hc + 10L, x 

(M = 1+ 47r W (-2L MS + 6L MX ), 

CM = 1+ ^ (-^-3^ + 4^) , (4) 

where L^ x = ln(yU, 2 /Mj) and for simplicity we keep from the list of arguments of the 
coefficients ( a% only the decoupling scale. 

The Missing Doublet Model [281129] is designed to avoid unnatural doublet-triplet splitting 
in the field H that is present in the minimal model. This is achieved at the cost of 
introducing additional Higgs fields and in the large SU(5) representations 50 and 50 
that do not contain any isodoublets and thus only couple to the colour triplets in H and 
H . To break SU(5) another Higgs field E in the 75 representation is used instead of the 
24 as in the minimal model. The superpotential reads 

W = M!Tr(S 2 ) + AiTr(S 3 ) + A 2 HS0 + A 2 HS0 + M 2 00 

+ v'2v;;M',o ; // + ^YzViVjH . (5) 
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After £ develops a vacuum expectation value the spectrum of the theory can be 
parametrized by five mass parameters [27] M x , M Hc , Mh c , , and M 2 . The last is as- 
sumed to be of 0(Mp{) so that the representations 50 and 50 do not contribute to the 
running above the unification scale and below Mj>\. Otherwise, due to large group factors 
of these representations, the perturbativity requirement cannot be fulfilled. In this case 
the decoupling constants read [26JU 



a sv( - 5 Uu) (22 64 
( ai {(i) = 1 + f -- L^ Hc - - L^ Hc , + 10 L^ x - 20 L m£ + 10 In — 

CaM = 1+ An KN f-22L MS + 6L^ + 61n- 



a sv ^(u) ( 64 \ 

C Q3 0*) = 1 + ^ \ -Wn c - L, Hcl - 23 + 4 L, x + 4 In _j . (6) 

As we will see later, the presence of the large representation 75 in this model leads to 
huge theoretical uncertainties due to the variation of the unphysical scale /iqut- 



2.2 Running and decoupling 

It is well known that gauge coupling unification is highly sensitive to the super-heavy 
mass spectrum [30]. This property allows us to probe unification through precision mea- 
surements of low-energy parameters like the gauge couplings at the electroweak scale 
or the supersymmetric mass spectrum. A simple algebraic exercise taking into account 
the naive step-function approximation [31] based on one-loop RGEs provides analytical 
formulas for the determination of the GUT spectrum as a function of the three gauge 
couplings measured at the Z-boson mass scale. An estimate of three-loop as compared to 
two- loop running has been obtained in Ref. [32] • Note, however, that a consistent treat- 
ment requires the implementation of ra-loop RGEs and (n — l)-loop threshold corrections. 
We have adopted this approach for n — 1, 2 and 3, whenever the required theoretical 
input was available, and have solved numerically the system of differential equations. 

Crucial input for our analysis constitute the precise values of the gauge couplings at the 
electroweak scale. They are obtained from the weak mixing angle in the MS scheme [15], 
the QED coupling constant at zero momentum transfer and its hadronic [32] contribu- 
tion in order to obtain its counterpart at the Z-boson scale, and the strong coupling 

2 The occurrence of the last term in the round brackets of each equation is due to the use of relations 
between the super-heavy masses. 
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constant [TO] HI The corresponding central values and uncertainties read 



sin 2 9 MS 

a 

a s {M z ) 



0.23119 ±0.00014 
1/137.036, 
0.02761 ±0.00015. 
0.1184 ±0.0020. 



(7) 



Whereas sin 2 9 MS and a s (Mz) are already defined in the MS scheme, Aa^j constitute 
corrections to the on-shell value of a. In order to obtain the corresponding MS result 

(5) 



we add the leptonic [34] and top quark [35] contribution, Aa| e p = 314.97686 x 10 and 



-0.70 ± 0.05) x 10 , and apply the transition formula to the MS scheme [15 



2 7 



1_7 M| 

e _ 4 ln Mi; 



This leads to 

a m {M z ) 



o 



1 - Act 



(5) 
lep 



Act 



(5) 
had 



Act 



(5) 
top 



0.0072 



0.0072 



127.960 ±0.021 



(9) 



In the quantities sin 2 6 MS , a MS (M z ) and a s (M z ) the top quark is still (partly) decoupled. 
Thus, in a next step we compute the six-flavour SM quantities using the relations [1~5"|I3"6] 



a 



(6), MS 



MS 



sin 2@(6),MS 



Oi 



sin 2 6 




where all couplings are evaluated at the scale jx = M z . We obtain0 



„(6),MS 



(M 2 



sm 2 0^ m (M z ] 



a 



(6) 



[M 2 



1/(128.129 ±0.021] 
0.23138 ±0.00014, 
0.1173 ±0.0020. 



3 We adopt the central value from Ref. [16], however, use as our default choice for the uncertainty 
0.0020 instead of 0.0007. 

4 Since we aim for gauge couplings at the electroweak scale with highest possible precision we use 
four-loop running and three-loop decoupling as implemented in RunDec [37j in order to obtain ai 6 ' from 
a s (Mz) = ai ,{Mz)- At such high order in perturbation theory there is practically no dependence on 
the decoupling scale. 
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These quantities are related to the three gauge couplings via the equations 



5 «( 6 )- MS 
3 cos 2 6(6)^ 



(6), MS 



sin 2 e(6)- MS 




(12) 



which holds for any renormalization scale fi. 

To the accuracy we are aiming at one has to worry about supersymmetric effects influ- 
encing the extraction of the couplings in Eq. ( ITT]) from experimental data. Due to the 
presence of the weak gauge bosons in the loop corrections the weak mixing angle receives 
the numerically largest contributions whereas the influence of supersymmetric particles 
on the electromagnetic and strong coupling can be neglected. 

The procedure to incorporate the supersymmetric effects on the numerical value of 
gin 2 e(6),MS( Mz ) ig ag follows: In a first step we transfer sin 2 9 (6) ' M5 (M Z ) from Eq. (JIT} to 
the DR scheme [38J and apply afterwards the supersymmetric one-loop corrections eval- 
uated in Ref. [39] relating the weak mixing angle in the SM to the one in MSSM. In a 
next step we decouple the supersymmetric particles [40] and finally go back to the MS 
scheme. As a result we obtain sin 2 Q^' MS (Mz) at the scale /i = Mz including virtual 
MSSM contributions. Note that by construction these corrections are suppressed by the 
square of the supersymmetric mass scale. We anticipate that for typical supersymmetric 
benchmark scenarios the influence of supersymmetric corrections to sin 2 Q^' MS (Mz) can 
lead to shifts in Mh c which are of the order of 10%. For the mSUGRA parameters in Eq. f fT6|) 
the shift of sin 2 6 (6) ' M5 (M Z ) amounts to about 1.4 ■ 10" 5 . 

In addition to the input values for the gauge coupling constants, we also need the W- 
and Z-boson pole masses My/ and M^, the top-quark and tau-lepton pole masses M t and 
M T and the running bottom-quark mass m^ s . For the convenience of the reader we also 
specify their numerical values [15ll4"Tll4*2] 



The corresponding uncertainties are not important for our analysis. 

In the following we describe in detail the individual steps of the running analysis needed 
for the energy evolution of the gauge couplings from the scale \i = Mz to the unification 




M w 
M z 
M t 
M T 



80.398 GeV, 
91.1876 GeV, 
173.1 GeV, 
1.77684 GeV, 
4.163 GeV. 



(13) 



scale. 
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1. Running within the SM from fi = Mz to the SUSY scale hsvsy- 

Starting from Eqs. (TTTT) and (fT2l (and adding supersymmetric effects as discussed 
above) we use the three-loop beta function of QCD [431144] and the two-loop RGEs 
in the electroweak sector [4"5T - T4"T] in order to obtain the values of the gauge couplings 
at /^susy ~ 1 TeV. We take into account the tau, bottom and top Yukawa couplings 
and thus solve a coupled system of six differential equations. Since the quartic 
Higgs coupling A enters the equations of the Yukawa couplings starting from two- 
loop order only we neglect its contribution. At this point we want to stress that 
A*susy is not fixed but kept as a free parameter in our setup. 

2. SUSY threshold corrections. 

For energies of about //susy ~ 1 TeV the SUSY particles become active and the 
proper matching between the SM and the MSSM has to be performed. We decouple 
all heavy non-SM particles simultaneously at the scale /isusy using the one-loop 
relations for a\ and a 2 and the Yukawa couplings from Refs. [10] and [38], respec- 
tively. The SUSY-QCD decoupling effects for and are known to two-loop 
order and have been computed in Refs. [TM19]. The simultaneous decoupling might 
be problematic in case there is a huge splitting among the SUSY masses. In that 
case a step-by-step decoupling would be preferable (see, e.g., Ref. [33]), however, 
a two-loop calculation in that framework is still missing. Furthermore, the mass 
splitting in almost all benchmark scenarios currently discussed in the literature is 
rather small. 

As pointed out before, a fully consistent approach would require two- loop threshold 
corrections not only in the strong but also in the electroweak sector. They are not 
yet available, however, we also expect that their numerical impact is relatively small. 
Furthermore, consistency with the RG running would require mixed QCD- Yukawa 
corrections for a% which are not yet available. Since the effect of these kind of 
corrections on the (three-loop) running is numerically small we expect that also the 
impact on the decoupling is small. 

In our numerical analysis we generate the SUSY mass spectrum with the help of the 
program SOFTSUSY [50] and study the various SPS (Snowmass Points and Slopes) 
scenarios [5TH52] . 

At this stage also the change of renormalization scheme from MS to DR has to be 
taken into account. We employed the one-loop conversion relations [38j|53] for all 
parameters except «3 and mj where two-loop relations [531154] have been used in 
order to be consistent with the decoupling at the SUSY scale. 

3. Running within the MSSM from /isusy to the high-energy scale /Uqut- 

We use the three-loop RGEs of the MSSM [I7J[20] to evolve the gauge and Yukawa 
couplings from /^susy to some very high scale of the order of 10 16 GeV, that we 
denote by fiQVT, where we expect that SUSY-GUT particles become active. 

4. SUSY-GUT threshold effects. 



S 



At the energy scale /xgut, threshold corrections induced by the non-degenerate 
SUSY-GUT spectrum have to be taken into account. Currently, only the one- 
loop corrections are available [22TI24"] . which we include in our setup, although for 
consistency two-loop corrections would be necessary. 

A suitable linear combination of the three one-loop equations for ccj in ([3]) and (j4]) 
leads to the following two relations 



4tt 



1 



4?r 5 



«iO) 
1 



a 2 (n) 
1 



- 2 



1 



12 



L 



)iH c i 



oli{h) a 2 {n) a 3 (n) 



-24 ( L, x + -L^ ) , (14) 



where a su ^ 5 ^ has been eliminated. These equations allow for the prediction of the 
coloured triplet Higgs boson mass Mh c from the knowledge of the MSSM gauge 
couplings at the energy scale \i = /igut- It is furthermore common to define a 
new mass parameter Mq = -^/M^Ms, the so-called grand unified mass scale, that 
can also be determined from the knowledge of the MSSM gauge couplings at /xgut- 
These observations makes it quite easy to test the minimal SUSY SU(5) model 
once the required experimental data are available in combination with a high-order 
analysis. 

For the Missing Doublet Model the above relations read 



Air 



1 



1 



4tt 5 



1 



-3 



1 



a 3 (/x) 
1 



Oi 2 (fx) a 3 (/i) 



12 , N , 64 

— (L„ Hc + L^) + 12 In 



-24 ( L^x + 



12 In 



3125 ' 
262144 
1953125 ' 

(15) 



5. Running from /igut to the Planck scale M P1 . 

The last sequence of our approach consists in the running within the SUSY-SU(5) 
model. We implemented the three- loop RGEs for the gauge [55], and the one-loop 
formulas for the Yukawa and Higgs self couplings [56] and impose that they can be 
described within perturbation theory up to the Planck scale. 



In Tab. [T] we summarize for the individual steps of the running-decoupling procedure to 
which loop order the perturbative corrections are currently available and implemented in 
our setup. The numbers in parenthesis indicate the loop order which would be necessary 
in order to perform a fully consistent analysis with three-loop running and two-loop 
decoupling. The phenomenologically most important ingredient which is still lacking 
are the two-loop decoupling relations at /^gut as we will show in the next Section. 
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run 

Mz — > ("SUSY 


dec 

(MSSM) 
OH -)• 0L\ 


run 

A^SUSY —> A^GUT 


dec run 

(MSSM) (GUT) > „ 


EW 
QCD 


2(3) 
3 


1(2) 
2 


CO CO 


1(2) 3* 



Table 1: Loop corrections available for the individual steps of the running-decoupling pro- 
cedure. The number is brackets indicate the loop-order needed for a consistent analysis. 
(* The running of the gauge and Yukawa coupling is performed to three- and one-loop 
accuracy, respectively.) 

3 Phenomenological constraints on the SUSY GUT 
mass spectrum 

The running-decoupling prescription described in Section 12.21 introduces two decoupling 
scales, /UsusY and //gut, for the supersymmetric and GUT particles, respectively. These 
scales are not determined by theory. However, on general grounds one expects that the 
GUT parameters become insensitive to the precise choice of these scales when going to 
higher orders. In Fig. Ufa) this is demonstrated for Mh c considering the dependence on 
Hsvsy where the latter is varied between 100 GeV and 10 TeV. For illustration we adopt 
the mSUGRA scenario for the SUSY breaking mechanism with 

m = mi/2 = -A = 1000 GeV , 
tan p — 3 , 

fi > 0, (16) 

and generate with the help of S0FTSUSY [50] the supersymmetric mass spectrum. This 
results in squark masses which are of the order to 2 TeV and thus in the upper range of 
what can be measured at the CERN LHC. The dotted, dashed and solid lines in Fig. [TJ^a) 
correspond to the one-, two- and three-loop running analysis, respectively, where the de- 
coupling is performed at one order lower as required by consistency. One finds a quite 
sizeable dependence of Mh c at one-loop order varying by almost two orders of magnitude. 
A significant reduction is observed after inclusion of the two-loop effects leading to a vari- 
ation of Mh c by only a factor two to three. Finally, after incorporating three-loop running 
and two-loop matching corrections the variation of Mh c on /^susy is about 5 ■ 10 14 GeV 
in the considered range for the matching scale. It is furthermore remarkable that for 
/U-susy around 1000 GeV, which is close to the average of the supersymmetric masses, the 
two-loop corrections show a maximum and the three-loop corrections are practically zero. 
In particular, they are significantly smaller than for /^susy = Mz which has often been 
used as decoupling scale for the supersymmetric particles. 

In Fig. [T](b) the dependence of M Hc on /igut is studied. A consistent three-loop analysis 
can not be performed since the two-loop GUT matching relation is not yet available. 
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Figure 1: Dependence of M Hc on //susy (a) and //gut (b) including successively higher 
orders. Dotted, dashed and solid lines correspond to the one-, two- and three- loop running 
analysis, respectively. 



Nevertheless it is tempting to combine the three-loop running with the one-loop matching 
effects which is represented by the solid line. One observes that Mh c varies by about 
1.5 • 10 15 GeV. This is of the same order of magnitude as the three- loop effect at the 
SUSY scale if the matching is performed at Mz (see Fig. QJa)). 

In the following we discuss the dependence of Mh c and Mq on various parameters entering 
our analysis. We start with varying the supersymmetric mass spectrum by considering 
different SPS scenarios [5"T||52] and use Eq. (fH|) in order to extract in each case both M Hc 
and Mq. The decoupling scales are fixed to /^susy = 1000 GeV and /igut = 10 16 GeV, 
respectively, which ensures, according to the previous discussion, that the three-loop effect 
is rather small. In Fig. [2] the results are shown in the Mh c — Mq plane where the lines 
indicate the slopes in case of SPSla, SPS2, SPS3, SPS7, SPS8 and SPS9. Most scenarios 
lead to M Hc masses between 0.25 • 10 14 GeV and about 1 • 10 15 GeV: SPS9 (anomaly- 
mediated SUSY breaking) gives the smallest and SPS2 the largest value of Mh c . 

A further illustration of the dependence of the GUT parameters on the SUSY spectrum 
can be found in Fig. [3] where we adopt the parameters of Eq. (fT6l) and vary mi/2- In 
this way we alter the SUSY spectrum which enters the prediction of Mh c and Mq via the 
decoupling procedure at \i = yUsusY- The solid and dashed lines correspond to M Hc and 
Mq, respectively, which show a substantial variation. On the other hand, m , tan/3 and 
A have only a minor influence on the GUT masses and thus we refrain from explicitly 
showing the dependence. 

In the following we study the effects of the uncertainties of the input values a.i (c.f. 
Eq. (Tl2l ) on Mh c and Mq. We fix the SUSY spectrum as before (see Eq. (fl6l) ) and set in 
a first step /isusY = Mz which has often been common practice in similar analyses (see, 
e.g., Ref. [9]). Taking into account correlated errors and performing a \ 2 analysis leads 
to ellipses in the M Hc — M G plane. Let us mention that we can reproduce the results 
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M H (10 15 GeV) 



Figure 2: Dependence of Mh c and Mq on the various SPS scenarios. 

of Ref. [9] after adopting their parameters and restricting ourselves to the perturbative 
input used in that publication. 

In Fig. ID^a) we show our results for the two- (dashed lines) and three-loop (continuous 
lines) analyses. In each case the two concentric ellipses correspond to 68% and 90% 
confidence level, respectively, where only parametric uncertainties from Eq. fTTTl) have been 
taken into account. A significant shift to higher masses of about an order of magnitude is 
observed for Mh c ; Mq increases by about 2 ■ 10 15 GeV. This demonstrates the importance 
of the two-loop matching and three-loop running corrections. As has been discussed in 
the context of Fig. CD they are essential in order to remove the dependence on /isusy- In 
fact, adopting in Fig. [4] //susy = 1 TeV would lead to ellipses for the two- and three- 
loop analysis which were almost on top of each other and which would coincide with the 
three-loop ellipses (solid lines) in Fig. HI 

As expected, the uncertainty of a s induces the largest contributions to the uncertainties 
of Mh c and Mq. In particular, it essentially determines the semimajor axis of the ellipses. 
Thus it is tempting to assign a more optimistic uncertainty of 5a s = 0.0010 and redo the 
previous analysis. The results, shown in Fig. 11(b), underline even more the importance 
of the three-loop analysis. Whereas for 5a s = 0.0020 the parametric uncertainty (i.e. the 
size of the ellipses) and the shift due to higher perturbative corrections are of the same 
order of magnitude, in Fig. H(b) the latter is about twice as big as the former. Let us, 
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Figure 3: Dependence of Mh c (solid) and M G (dashed) on m^. 



however, stress once again that choosing /isusY close to the supersymmetric mass scales 
leads to small three-loop effects since the two-loop ellipses are essentially shifted on top 
of the three-loop ones. 

Recently there have been a few extractions of a s based on higher order perturbative 
corrections with uncertainties slightly above 1%, which, however, obtain central values 
for a s close to 0.113 (see, e.g., Ref. [5?]). Since these results are significantly lower than 
the value given in Eq. (J7J) it is interesting to show in Fig. |4Ja) also the corresponding 68% 
and 90% confidence level ellipses (for the two- and three- loop analysises) as dotted lines 
adopting a s [Mz) — 0.1135 ± 0.0014 [58]. One observes a big shift in the GUT masses, in 
the case of M Hc the central value is about one order of magnitude lower than for the a s 
value of Eq. ©. 

In Fig. [5] we visualize the running (and decoupling) of the gauge couplings where the 
parameters of Eq. ffTB"]) together with /^susy = 1000 GeV and /iout = 10 16 GeV have been 
adopted. In addition we have chosen = 1-10 15 GeV which leads to Mh c = 1.7-10 15 GeV 
and Mx = 4.6- 10 16 GeV. One can clearly see the discontinuities at the matching scales and 
the change of the slopes when passing them. In panel (b) the region around /i = 10 16 GeV 
is enlarged which allows for a closer look at the unification region. The bands indicate la 
uncertainties of ccj at the electroweak scale (cf. Eq. (ITT]) ). In panel (b) we furthermore 
perform the decoupling of the super- heavy masses for two different values of /igut- One 
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13 14 15 16 

log 10 (M Hc /GeV) 

(b) 

Figure 4: Ellipses in the Mh c — Mq plane obtained from the uncertainties of the gauge 
couplings at the electroweak scale. In (a) the input parameters of Eq. (J7J) have been used 
whereas in (b) 5a s = 0.0010 has been chosen. Dashed and solid lines correspond to the 
two- and three-loop analysis, respectively. The dotted lines in (a) have been obtained for 
a s (Mz) = 0.1135 ±0.0014 where the lower (upper) ellipse corresponds to the two- (three-) 
loop analysis. 
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Figure 5: Running of the gauge couplings from the electroweak to the Planck scale. The 
discontinuity for \i = /isusy and /i = /xqut are clearly visible. In panel (b) an enlargement 
of (a) for the region around fi = //gut is shown where for the decoupling the two values 
A*gut = 10 15 - 5 GeV « 3.2 ■ 10 15 GeV and /i GU t = 10 16 GeV have been chosen. 

observes quite different threshold corrections leading to a nice agreement of a su ^ 5 ^ above 
10 16 GeV. Fig. [5] stresses again that the uncertainty of a s is the most important one for 
the constraints that one can set on GUT models from low-energy data. Furthermore, it 
illustrates the size of the GUT threshold corrections and emphasizes the importance of 
the two-loop corrections for the corresponding decoupling constants. 

Finally we discuss the phenomenological consequences of our analysis in a "top-down" 
approach where we specify the parameters at the high scale and examine the effect on 
the gauge couplings for \x = Mz- For minimal SUSY SU(5) we choose the following 
parameters 



which guarantee the agreement of a (6) - MS (M z ), sin 2 6( 6 )- MS (M z ) and af )M (M z ) with 
their experimental counterparts for //gut = 10 16 GeV. In Fig.[6ta) we fix /xsusy = 500 GeV 
and the SUSY spectrum according to the SPSla benchmark scenario which allows for a 
comparison with the MDM (cf. Fig. EJ^b)). We vary /^gut by three orders of magnitude 
where the dotted, dashed and solid line correspond to the one-, two- and three-loop 
analysis of the described procedure, respectively, and the bands reflect the uncertainties 
as given in Eq. (II ip . For this parameter choice we observe, as expected from the above 
discussions, a sizeable two-loop effect but only a mild change after including the three- 
loop corrections. Furthermore, the overall variation is small for all three gauge couplings. 
Note that the one-loop curve for a s is independent of //gut since the one-loop coefficients 
of the (3 functions in SUSY QCD and minimal SUSY SU(5) are identical. 



M x 



3.67 • 10 14 GeV 
2- 10 16 GeV, 
1.58 ■ 10 16 GeV 



a 



su ( 5 )(10 17 GeV) 



0.03986 



(17) 
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In Fig. EJb) we perform the same 
500 GeV and the parameters^] 



analysis within the MDM adopting SPSlaJfl /isusy 



a 



su ( 5 )(10 17 GeV) 



M x 



6-10 18 GeV, 
1 • 10 16 GeV, 

2- 10 15 GeV, 

3- 10 16 GeV, 
0.1504. 



(18) 



One observes a much stronger variation of the predicted values of a^, in particular in the 
case of the strong coupling which varies between 0.06 and 0.13 in the considered range of 
A*gut- This demonstrates the importance of the two- loop threshold corrections which are 
expected to significantly reduce the dependence on //gut- 

At this point a discussion about the additional constraint on the Higgs triplet mass 
Mh c that can be derived from the non-observation of the proton decay is in order. The 
latest upper bound on the proton decay rate for the channel p —> K + u [TJ] is r exp = 
4.35 x 10~ 34 /y. In order to translate it into a lower bound for the the Higgs triplet mass, 
one needs an additional assumption about the Yukawa couplings that enter the expression 
of the decay rate F(p — > K + u). As pointed out in Ref. [11] this is because down quark 
and lepton Yukawa couplings fail to unify within the minimal renormalizable SUSY SU(5) 
model and so a completely consistent treatment is not possible. Therefore one could either 
choose^ (i) Y ql = Y ud = Y d or (ii) Y q i = Y ud = Y e , which leads to completely different 
phenomenological consequences. We are aware that both cases are equally justified once 
higher dimensional operators are included. Since these operators further weaken the 
bounds presented below, we refrain to include them into the analysis in this paragraph. 
For the case (i) and sparticle masses around 1 TeV the lower bound for the Higgs triplet 
mass can be read off from Fig. 2 of Ref. [TT] and amounts to Mh c > 1 .05 x 10 17 GeV whereas 
for the second choice it becomes M# c > 5.25 x 10 15 GeV. From our phenomenological 
analysis presented above it turns out that within the minimal SUSY SU(5) model the 
upper bound for Mh c is of about 10 16 GeV. So, the substantial increase of about one order 
of magnitude for the upper bound on M Hc induced by the three-loop order running analysis 
attenuates the tension between the theoretical predictions made under the assumption (i) 
and the experimental data. The choice (ii) for the Yukawa couplings clearly shows that 
the minimal SUSY SU(5) model cannot be ruled out by the current experimental data on 
proton decay rates. More experimental information about the SUSY mass spectrum and 
proton decay rates is required in order to be able to draw a firm conclusion. 

Adopting the parameters from Eq. (|16j) leads to a non-perturbative values of the gauge coupling at 
the Planck scale. 

6 The relatively large value of a su ( 5 )(10 17 GeV) is required due to the large Casimir constants entering 
the /3 function above the GUT scale. 

r Y q i is the Yukawa coupling of the quark and lepton doublet to the Higgs colour triplet. Y u d is the 
corresponding coupling for the up and down quark singlet. 
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(b) 

Figure 6: Gauge couplings at the weak scale as obtained by a top-down approach within 
minimal SUSY SU(5) (a) and the MDM (b) as a function of /^gut (see text). In (b) the 
one-loop curves are only shown up to yUcuT ~ 4 ■ 10 14 GeV since beyond this scale 
becomes quite large. 
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On the other hand, it is commonly recognized that SUSY SU(5) in its minimal version 
cannot provide the underlying model for a consistent GUT because of the failure of Yukawa 
unification. As mentioned in the Introduction, one of the most appealing solutions of this 
problem is the inclusion of Mpi scale effects using the effective theory approach, i.e., 
including higher dimensional operators. Explicitly, this translates into [T3] 



where M H denotes the mass of the Higgs triplet in the minimal renormalizable model. 
In this case, a suitable choice of the Yukawa matrices can significantly reduce the proton 
decay amplitudes, so that values for M Hc of the order C(10 15 GeV) even for larger values 
of tan /3 ~ 15 are allowed. So the prediction of the non- renormalizable version of the 
SUSY SU(5) model is well under the experimental upper bound. However, the inclusion 
of new undetermined parameters together with the higher dimensional operators makes 
the tests of such models much more complicated. 

At this point we try to compare with the findings of the analysis [59] based on the non- 
renormalizable version of minimal SU(5) where the bound Mh c > 3.7 x 10 17 GeV is 
given. This is the most stringent bounqj that can be derived from current proton decay 
experimental data. In order to be consistent with our upper bound on Aff 5 (cf. Fig. H} 
of about 10 16 GeV one requires M(i^/M(s ^ > 4.2 which is a quite moderate value. 

4 Conclusions 

Already in the early days of the Standard Model there have been studies of theories which 
predict the unification of the coupling constants at high energies. Since it is not possible 
to reach such energies in collider experiments it is necessary to establish relations between 
the couplings at the electroweak scale, where precise measurements are available, and the 
corresponding quantities at the unification scale. Next to the beta functions covering the 
running also threshold corrections at the SUSY and GUT scale constitute crucial input. 

We have considered the so-called minimal supersymmetric SU(5) GUT theory and have 
studied the gauge coupling unification applying one-, two- and three- loop running. The 
main effect of the three-loop running (accompanied by two-loop decoupling relations) is 
the stabilization w.r.t. the variation of the decoupling scales. In general sizable three- 
loop effects are observed if the decoupling scale for the supersymmetric particles is not 
chosen in the vicinity of the masses of the supersymmetric particles. In particular, for 
A*susy = Mz, which is the canonical choice often adopted in the literature, one observes 
an increase of the coloured Higgs triplet mass by about an order of magnitude. Already in 
previous studies it has been shown that the non-renormalizable version of minimal SUSY 

8 See Ref. [59] for the exact definition of this bound. 
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(19) 
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SU(5) cannot be excluded by the experimental bound on the proton decay rate [T0l[TTl[T3] . 
Our results attenuate the tension between the theoretical predictions of Mh c and the 
experimental results even more. 

Our analysis includes all the state-of-the-art theoretical input. It could be improved 
by including the two-loop GUT threshold corrections which are not yet available. This 
induces an uncertainty of about 0.3 on log 10 (M Hc /GeV) which constitutes the major 
theory-uncertainty. We observe that the effects of the two-loop GUT threshold corrections 
become particularly important in the Missing Doublet Model where also larger effects can 
be observed. As far as the parametric uncertainty is concerned we mention the dependence 
on the supersymmetric mass spectrum and the uncertainty on a s . In both shift 
of Mh c of a few times 10 15 GeV is observed. 

The most popular extensions of the minimal SUSY SU(5) model, the Missing Doublet 
Model and the non-renormalizable version of the SUSY SU(5) model comprising Planck- 
scale operators, can not be excluded using only the currently available theoretical and 
experimental data. They either contain additional free parameters as compared to the 
minimal model or are affected by large theoretical uncertainties, so that no firm conclusion 
can be drawn for them. 
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